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Abstract 



The Lie-group approach to the perturbative renormalization group 
(RG) method is developed to obtain an asymptotic solutions of both 
autonomous and non-autonomous ordinary differential equations. Re- 
duction of some partial differetial equations to typical RG equations is 
also achieved by this approach and a simple recipe for providing RG 
equations is presented. 
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1 Introduction 



A novel method based on the perturbative renormalization group (RG) the- 
ory has been developed as an asymptotic singular perturbation technique 
and has been succesively applied to ordinary differential equations (ODE) |2 
and some partial differential equations (PDE) @] [0 f| • An approach from 
a geometrical point of view was also proposed, which is called an envelope 
method 0. 

The renormalization group method removes secular or divergent terms from a 
naive perturbation series by renormalizing integral constants i.e. a renormal- 
ization transformation and reaches a slow- motion and/or large-scale equation 
as a renormalization group (RG) equation. In this reduction of an asymptotic 
equation, it is not necessary to introduce slow or large scales a priori in con- 
trast to the other methods such as the reductive perturbation method |9j and 
various multiple-scale methods [10]. Instead, a renormalization transformation 
determines an asymptotic equation as a RG equation. Here, it is crucial to 
extract the structure of secular terms from a perturbation series in order to 
obtain a renormalization transformation. For ODE, the structure of secular 
terms are often relatively simple and, as shown in this paper, it is possible 
to calculate explicit secular solutions up to arbitrary order in principle and 
obtain an asymptotic form of the renormalization transformation of all order 
in some examples . However, for PDE, there are many divergent solutions in 
a perturbed equation and it is ambiguous which divergent solutions should be 
renormalized away M. Although an attempt to overcome this difficulty was 
presented ||, it is far from satisfactory. 

The group structure of the renormalization transformation has not yet been 
necessarily clear in the previous applications of the RG method and it is desire- 
ble to make the renormalization procedure more transparent for its extended 
applications. 

In this paper, we construct a representation of the Lie group from a renor- 
malization transformation in several examples and identify the procedure of 
renormalization with one to derive an asymptotic expression to a generator 
of the Lie group. Although the view of the Lie group has been introduced to 
the RG method in systems with a translational symmetry with respect to a 
independent variable [11|@], there are few examples based on the Lie group 



theory and only the leading or low order RG equations are derived || . For the 
purpose of obtaining higher-order RG equations and applying the RG method 
to wider systems such as non-autonomous systems without a translational 
symmetry and PDE, much more investigations are still necessary to refine the 
Lie-group approach. 

In Section 2, we present some clear examples, in which an asymptotic ex- 
pression to a generator of the Lie group (a renormalization transformation) 
is obtained up to arbitrary order, in principle, by naive pertubation calcula- 



2 



tions. In the simplest example in Section 2. A , an asymptotic expression of 
a generator is shown to converge and we recover an exact solution. In Sec- 
tion 2.B , we extend the Lie-group approach to a non-autonomous system of 
ODE by introducing a simple shift operation. Through these simple exam- 
ples, we present a recipe for obtaining an asymptotic representation to the 
Lie group underlying the original system. In the remaining parts of Section 
2, we give some physically interesting examples. In an example in Section 
2.C, a renormalization transformation has a translational symmetry though 
an original ODE does not have one. In Section 2.D, we present a peculiar 
example in the Einstain's gravitational equation where an original ODE and 
its leading order RG equation are autonomous but a higher order RG equation 
becomes non-autonomous. Even in this peculiar case, our Lie-group approach 
works well and gives an asymptotic solution describing an expanding space. 
In Section 2.E, it is shown that secular solutions to perturbed equations for 
a weakly nonlinear oscillator can be removed asymptotically up to all order. 
We also derive adiabatic equations for parameters of the K-dV soliton under 
perturbations in Section 2.F. 

In Section 3, an extension of the RG method in terms of the Lie group is 
presented for some nonlinear PDE. We derive two typical soliton equations 
such as the nonlinear Schrodinger equation and the Kadomtsev-Pitviashvili 
equation with higher order corrections as RG equations. As the final exam- 
ple, phase equations for interface and oscillating solutions are derived from a 
general react ion- diffusion system. 

In Section 4, we give a simple recipe for providing RG equations as a summary 
of this paper. 



2 Ordinary Differetial Equations 



2. A Linear Autonomous System 

As an example in which a generator to the Lie group of a renormalization 
transformation is constructed exactly, we discuss a boundary layer type prob- 



where u(0) = 0. Since a small parameter e multiplying the highest derivative 
causes a severe singularity at a boundary x = 0, we introduce a transformation 



lem irrn . 
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et. Then, ( |2.1| ) reads 



( d 2 d \ tnns 

{w + T t + j u = < (2 ' 2) 

which is called the inner equation in the boundary layer problem. Expanding 

u as 

u = A + Be~ l + em + e 2 u 2 H , (2.3) 

we have 

/ d 2 d \ 

L(t)u n =^- + -)u n = -u n „ u (2.4) 

where uq = A + Be~ l and A, B are arbitrary constants. Secular solutions can 
be written as 

u n = AP n (t) + BP n {-t)e-\ (2.5) 
where Pq = 1 while P n {t) (n > 1) is a polynomial solution of 

L(t)P n = — P n -i, 

P 1 = -t, P 2 = t 2 /2-t, P 3 = -t 3 /6 + t 2 - 2t. 

Here, we observe that a secular term P n (t) is a polynomial of degree n. In order 
to eliminate these secular terms by renormalizing integral constants A, B, the 
following naive renormalization transformations A — > A(t) and B — > Bit) are 
introduced as 

A{t) = A(l + eP 1 (t) + e 2 P 2 (t) + •••), (2.6) 
Bit) = B{1 + eP 1 {-t) + e 2 P 2 {-t) + ■■■). (2.7) 

Note that this definition of renormalization transformations is different from 
one in 0. In order to obtain a representation to the Lie group from ( |2.6| ). 
we take advantage of a translational symmetry in ( f2.2p . It is easy to see that 
A(t),B{t) also enjoy a translational symmetry ,that is 

L(t)A + eA = 0, L(-t)B + eB = 0. (2.8) 

The renormalization transformation (|2.6|) is interpreted as a Taylor series for 

the solution A of ( |2.8| ) around t = 0, while the origin of expansion can be 
shifted to arbitrary positions by virtue of a translational symmetry. Thus, we 
obtain an asymptotic representation to the Lie group G T immediately from 

& 

A(t + r) = A(t){l + eP 1 (r)+e 2 P 2 {r) + ---), 

= G T A{t). (2.9) 
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In terms of a generator rd t of the Lie group, A(t + r) is also expanded as 

A(t + r) = exp{Tdt)A(t), 

= (l + Td t + (r 2 /2)d t 2 + ---)A(t). (2.10) 

Equating coefficients of equal powers of r in djjD and (CTD , we have an 
asymptotic form of a generator of the Lie group called a RG equation 

HA 

^ = 9 T [G T i(t)] r=0 , 

= d T [eP 1 (r) + e 2 P 2 (r) + ---] T=0 A, 

= -e(l + e + 2e 2 + ---)A = eSA, (2.11) 

All other equations for higher order derivatives are guranteed to be consistent 
with (|2.11 ) by virtue of the theory of Lie group or arbitrainess of r. Since 
5 = -(1 + e + 2e 2 + • • •) converges to (-1 + ^1 -4e)/(2e) for e < 1/ 4 (ETTTfl 
gives an exact form of the generator. Applying a similar analysis to (|2.7), we 
obtain a RG equation for 5: 

^ = 9 r [ e P 1 (-r) + e 2 P 2 (-r) + ---] T=0J B, 

= e(l + e + 2e 2 + • • -)B = ~eSB. 

Thus,renormalizing all secular terms in the perturbation series, we recover an 
exact solution of ( |2.1| ) with the boundary condition u(0) = 0. 

u = 1(0) exp(5x) + 5(0) exp(-^ - 8x), (2.12) 
where i(0) + B(0) = 0. 



2.B Oscillator with time-dependent spring constant 



To extend the Lie-group approach described in the previous section to systems 
without a translational symmetry with respect to t, we consider an linear but 
nontrivial oscillator governed by 

d 2 u 

— + u = etu, (2.13) 
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which was analyzed in 0] to the first order in e. Here, we try to extend the 
RG analysis to arbitrary order in terms of the Lie group. Expanding u as 

u = Ae~ u + eu! + e 2 u 2 -\ I- c.c., (2.14) 

we have 

+ 1^J u n = tu n _ u (2.15) 

where c.c. stands for the complex conjugate of the preceding expression, n = 
1, 2, • • • and uo = Ae~ lt . Secular solutions are given by 

u n = AP n (t)e- u , 
L(t)P n (t) = (J^-2tj^JP n = tP n _ u (2.16) 

where n = 1, 2, • • • and Pq = 1, while an initial condition P n (0) = is set for 
n > 1. The important observation is that P n (t) is a polynomial of degree 2n 
which is given by, for example, 

P 1 = (tt 2 + t)/A, P 2 = y £ Pj iP, 

3=1 

p 4 = -1/32, p 3 = 5i/48, p 2 = 5/32, p x = 5z/32. 
Thus, we have 

u = A{1 + ePi(t) + e 2 P 2 (t) + ■ ■ .)e _il + c.c, (2.17) 

which yields a renormalization transformation A — > A(t) defined by the same 
expression as the previous example (|2.6| ), where secular terms P n are given by 
polynomial solutions of (|2.16| ). Since A(t) has no longer a translational sym- 
metry in this case, it is necessary to introduce the following general procedure 
in order to derive an asymptotic representation to the Lie group from ( |2.6| ). 
Let us shift t — > t + r in 



A(t + t) = A{1 + eP t {t + r) + e 2 P 2 {t + T ) + ■■■). (2.18) 

Noting that P n is a polynomial of degree 2n, we have the following expansion 
of A with respect to r and e. 



oo 

A(t + r) = A(l + ^e"P n (t)) + [r^e"P n , t 



n=l n=l 

oo oo 

+ (r 2 /2) £ e n P n , 2t + (r 3 /3!) £ e n P n , 3t 

n=l n=2 

oo 

+ (r 4 /4!) 5>™P n , 4t + ■ ■ (2.19) 

n=2 
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where P n ^t denotes the k-th derivative with respect to t. Replacing A in 
( P-19| ) with A by means of the renormalization transformation (|2.6| ), we reach 
an asymptotic representation to the Lie group : 



A(t + r) = A(t) + [rJ2^Pn,t + (r 2 /2)J2e n Pn,2t 

n=l n=l 

oo oo 

n=2 

oo 

+ H^(*)/(i + E^»(*))- 



n=2 



(2.20) 



n=l 



Due to loss of the translational symmetry, the right-hand side of ( j2.20 ) explic- 
itly depends on t and so the representation (|2.20|) should read 



A(t + r) = G T {A(t),T(t)}, T(t)=t. 



(2.21) 



From ( [2.20 ) and ( |2.21| ) , we can derive an asymptotic expression of a generator 
of the Lie group, that is, a RG equation: 



d A 

— = d T [G T {A(t),T(t)}]r=o 



d 
dt 



{hx(l + J2e n P n (t))}A(t), 



(2.22) 



n=l 



from which the following asymptotic solution A(t) is easily obtained by inte- 
grating a truncated expression of ( j2.22p . 



A(t) = A(0) exp{ePi(t) - e 2 (P 2 (t) + P 1 (t f) + ■••}, 



(2.23) 



where the argument of exp function in ( |2.23| ) is determined by integration of 
a truncated power series of 



dt 



(2.24) 



The solution ( |2.23|) is a genaralization of the previous result [0. Note that 
the RG equation ( 2.22j ) is formally integrated to recover the renormalization 
transformation (|2l 



Now, we summarize the above procedures for providing a RG equaion as 
a generator of a renormalization group (the Lie group). 
(1) Caluculate secular terms in a naive perturbation series. 
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(2) Find integral constants A E R n , which are renormalizable in order to 
remove secular terms, and construct a renormalization transformation A(t) = 
R(A, t) which is obtained as a power series of t. 

(3) Rewrite the renormalization transformation in the form of the Lie group 

A{t + t) = G T {A(t), T(t)} by means of a shift operation t — > t + r, where A 
and T(t) = t constitute a differentiate manifold on which the Lie-group G T 
acts. 

(4) A RG equation is obtained as a generator of the Lie group: 

d A 

— = d T [G T {A(t),T(t)}] T=0 . 

If a renormalization transformation A(t) = R(A, t) has a translational sym- 
metry with respect to t, a repreresentation of the Lie group does not depend 
on T(t) = t and takes the following simple form. 

A{t + T) =G T A(t) = R{A{t),T). (2.25) 

In the following subsections, we derive various RG equations following the 
above procedures. 



2.C Mathieu Equation 



As a second example of non-autonomous systems, we derive a RG equation 
near the transition curves that separate stable from unstable solutions of the 
Mathieu equation 

d 2 u 

+ u = - e {cj + 2cos(2t)K (2.26) 

at 2 

where ui = cu 1 + eu 2 + e 2 u 3 + • • ■. Expanding u just as in the previous example 

u = Ae~ u + c.c. + eui + e 2 u 2 H , (2.27) 

we have, in the leading order, 

a L^l + Ul = -{uiA + A*)e~ u - Ae~ m + ex., (2.28) 
at 1 

of which general solution is 

ui = -{uj\A + A*)P 1 (t)e~ i ' + A/8e- 3lt + c.c, (2.29) 
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where A* denotes the complex conjugate of A and P\{t) is a polynomial solu- 
tion of 

(£-*l)* s£W * =i - < 2 - 30 > 

where Pi(t) is uniquely determined for the initial condition -Pi(O) = 0, i.e. 
Pi{t) = it/2. Similar perturbation calculations yield higher order solutions: 

u 2 = -{{u 2 + l/8)P 1 + {l-oo 2 1 )P 2 }Ae- it 

+ {{u^A + A*)Qi + u} l A/8 2 }e~ m + (A/192)e- 5it + c.c, 

u 3 = [{(lu 1 /32-lu 3 )A + (3/6A)A*}P 1 + {2(lu 2 + 1/8) 

+ (l-u 1 2 )/2}Lu 1 AP 2 + (l-Lu 1 2 )(uj 1 A + A*)P 3 }e- it + c.c., (2.31) 

where, for simplicity, higher-harmonic terms in ( |2.31| ) are not listed explicitly 
and secular terms P 2 , P 3 , Q\ are polynomial solutions of 

LP 2 = Pi, LP 3 = P 2 *, 

Explicit forms of secular terms are given by 

P 2 = -t 2 /8 + it/8, P 3 = -it 3 /48, Qi = -it/16 - 3/64, (2.33) 

where Q\ is a non-resonant secular term in the third higher homonics, which 
comes from a resonant secular term P x in the fundamental harmonic. In order 
to remove secular terms in the coefficient of the fundamental harmonic (e~ rf ) 
to 0(e 3 ), we introduce a renormalization transformation 

A(t) = R(A, t) = A- e{u x A + A*)P X - e 2 {(u 2 + 1/8)P 1 + (1 - ujf)P 2 }A 
+ e 3 [{(^/32 - u 3 )A + (3/QA)A*}P 1 + {2(uj 2 + 1/8) 
+ (1-uj 1 2 )/2}uj 1 AP 2 + (1-uj 1 2 )(uj 1 A + A*)P 3 }. (2.34) 

Performing a shift operation t — > t + r on ( |2.34| ) and replacing A in the 
right-hand side of tfZM ) by A{t) + e(wii4 + A*)P x {t) + 0(e 2 ) , we obtain a 
representation of the renormalization group 

A(t + r) = G T A(t) = R(A(t),r). (2.35) 

This expression indicates that the renormalization group has a translational 
symmetry. In fact, ( |2.35|) is also directly derived from ( p.34| ) when we take 
into account of the fact that Fourier components of u obey a coupled but au- 
tonomous system of equations. This result shows an interesting example where 
the renormalization group has a translational symmetry even if the original 
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system does not have one. It is easy to see that non-resonant secular terms 
such as Qi in higher harmonics are automatically eliminated when resonant 
secular terms in the fundamental harmonic are renormalized by ( 12.34] ). Fur- 
ther discussions are presented about automatic elimination of non-resonant 
secular terms in higher harmonics in Section 2.E. By differentiating fl2.35|) 
with respect to r, we have a generator of the Lie group (|2.35| ) 

HA 

= -{ie/2){u!A + A*) - {te 2 /2)AA 

+ (ie 3 /2)[(uj 1 /32 - u 3 )A + (3/64)1* + (u 1 /2)AA], (2.36) 
A = uj 2 + 1/8 + (1 - w?)/4, 

from which 

^ = [ e 2 (l-^ 2 )/4-(eV)A 

+ e 4 {u 1 (u 1 /32 - u 3 ) -3/64 + WiA/2}]A 

which gives the transition curves that separate stable from unstable solutions: 
uj\ = ±1, co>2 = — l/8, c^3 = =Fl/64. This transition curves oji,u 2 and a 
special RG equation for u)\ — 1 were obtained to 0(e 2 ) by means of the RG 
method 0. The expression of a generator ( ]2.36| ) is a generalization of the 
previous results. 



2.D Large-scale expansion in Gravitational Equation 



A uniformly and isotropically expanding space with dust is discribed by the 
following Einstein's equation called the Friedman-Robertson- Walker equation 

2ad + a 2 = —(k - Aa 2 ), (2.37) 

where a is a time-dependent scale factor of space, a dot on a denotes the 
derivative with respect to time t, k(— ±1,0) is a sign of curvature and A is 
the cosmological constant. Although solutions of ( |2.37| ) with an initial condi- 
tion a(0) = are well-known and describe various kinds of expanding and/or 
contracting space, this example gives a peculiar asymptotic form of the repre- 
sentation of the Lie group for a large-scale solution as is shown later. That is, 
the higher order RG equation becomes non-autonomous though the original 
equation ( |2.37| ) and the leading order RG equation are autonomous. 
The leading order RG equation has been derived for a large-scale solution of 
( ]2.37| ) with A = by Y.Nambu and Y.Y.Yamaguchi [I2| and they shows that 
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a contraction stage of the Friedman space (k = 1) is qualitatively recovered 
by only the leading order RG equation. 

Here, we derive a higher order and non-autonomous RG equation by means 
of the present Lie-group approach for the purpose of quantitative comparison 
with the exact results. 

If A = 0, we can expand a around a large-scale solution ao = At 2 / 3 as 

a = x {A + ai{x)/A + a 2 {x)/A 3 + a 3 (x)/A 5 + ■••}, (2.38) 



where x = t 2//3 , A is an arbitrary constant and A > 1. Substituting ( |2.38|) 



into ( |2.37|) and solving the resulting equation order by order with respect to 



A 1 , we obtain the following secular terms 

aj = - Cj x j j = 1,2,3,...., (2.39) 

9k 3c 2 _ 23cic 2 _ 9c! + iQciCs 

Cl " 20' C2 "~' C3 "^T' C4_ 22 • 

Renormalizing the integral constant A so that secular terms ( |2.39| ) are re- 
moved, we have a renormalization transformation 

A(x) = A + ai {x)/A + a 2 (x)/A 3 + a 3 {x)/A 5 + ■■■. (2.40) 

Replacing x in (|2.40|) by x + £, we obtain a non-autonomous representation of 
the Lie group: 

A(x + t)=Gz{A(x),x} 1 
of which generator is explicitly given as 

dA 

— = <9 f G 5 (i(x),x)| f=0 

ci (cj - 2c 2 )x (2c? - 7cic 2 + 3c 3 )x 2 



A A 3 A 5 

(5cj - 24c 2 c 2 + 6c 2 + 16 Cl c 3 - 4c 4 )x 3 

+ ^ +•••. [ZAl) 

It is easy to see that this expression of the generator is asymptotic for x/A 2 = 
t/a(t) < 1, that is, the scale factor a is greater than the cosmic horizon. Re- 
sults of numerical integration of ( |2.41| ) for k — 1 are compared with the exact 
solution in Fig.(l). Solutions of the RG equation with higher-order corrections 
( |2.41| ) fits the exact solution closer in the region t/a(t) < 1 but they separate 



away from the exact solution for t/a(t) > 1 where the present asymptotic ex- 
pansion is not valid. It may be interesting that the leading order RG solution 
gives the best fit for t/a(t) > 1. This may come from fact that the leading 
order RG equation does not have the factor x/A 2 = t/a(t) , which should be 
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small in the present asymptotic expansion . 



When a value of A is near the critical value where the Friedman space 
begins to expand again, A should be scaled as 



A = X/Aq, (2.42) 

where Aq = A(0) and the critical value of A is given by A c = 9/4. Then, effects 
of this small A enter into (|2.41|) as the following correction terms of 0(A~ 5 ) 



and 0(A- 7 ) 

3X 3 Ax 2 + 2(ciA 3 - 2A 4 )x 3 , (2.43) 

where A 3 = \/(12A%), A 4 = 9Aci/(88Ag). Numerically integrating f^ill) with 
corrections (|2.43| ), we obtain asymptotic solutions near the critical value A c as 
shown in Fig. (2), which demonstrates again that asymptotic solutions agree 
well with an exact solution for t/a(t) < 1. 



12 




Figure 1: The scale factor a vs. t for A = the line is the cosmic horizon, 
i.e. a = t, the curve e is the exact solution, the curve 1 is a solution of the 
leading order RG equation, the curves 2 ,3 and 4 are solutions of the second 
and the third and the fourth order RG equations respectively. 
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Figure 2: The scale factor a vs. t for A = 1.1 A c the line is the cosmic 
horizon (a = t), the curve e is the exact solution, the curve 1 is a solution of 
the RG equation. 
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2.E Weakly Nonlinear Oscillator 



Let us apply the Lie-group approach to the motion of a weakly nonlinear 
oscillator. 

d 2 u 

—— + u = ea 3 u 3 + e 2 a 5 u 5 + e 3 a 7 u 7 H . (2.44) 

at 2 

Expanding u as 

u = Ae~ u + c.c. + eui + e 2 u 2 H , (2.45) 

we have 

+ 1 ) Ul= ^\A\ 2 Ae- lt + A 3 e~ 3li + c.c.]. (2.46) 
The general solution is given by 

Ul = a 3 [3\A\ 2 AP ie - tt + A 3 qi e~ 3it + c.c], (2.47) 



where q± = —1/8 and Pi(t) is a polynomial solution of (|2.30|) . For 0(e 2 ) , we 
have 

u 2 = {3a 2 3 {3P 2 + q 1 P 1 ) + 10a 5 P 1 }\A\ 4 Ae~ it 
+ {9a&i + (Qa 2 3qi + ba 2 )q x }\A\ 2 A 3 e~ m 
+ (3a 2 3 q x + a rj )q 2 A 5 e' 5it + c.c, 
LP 2 = P 1; fe = -1/24, L 3 Qi = Pi, (2.48) 

where L, L 3 are defined in (|2.30|) and ( |2.32|) . Explicit forms of secular terms 
are given in ( |2.33|) . Up to 0(e 2 ), the secular solution is given by 



u = [A + 3ea 3 \A\ 2 AP 1 + e 2 {3al(3P 2 + q 1 P 1 ) + 10a 5 P 1 }\A\ 4 A)e- it 
- (e/8){a 3 + e(9a^Pi + 21^/8 + ha 5 )\A\ 2 } A 3 e~ 3it + c.c. + ■ • • .(2.49) 



In (|2.49|) , coefficients of e %t are summed up to give an asymptotic expan- 
sion to the renormalized amplitude A . 

Ait) = A + 3ea 3 \A\ 2 AP l (t) + e 2 {3a 2 3 (3P 2 (t) + qi P x {t)) 

+ 10a 5 Pi(t)}|A| 4 A + O(e 3 ), (2.50) 

while coefficients of the third harmonic are summed up to give eA 3 [t) as 

eA 3 (t) = -(e/S)[a 3 + 9ea 2 3 \A\ 2 P l (t) 

+ e(21a 2 3 /8 + 5a 5 )\A\ 2 + 0(e 2 )}A 3 . (2.51) 
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Substituting an iterative expression of A in terms of A obtained from ( j2.50| ) 
into ( |2.51| ), we can eliminate a secular term in ( |2.51| ) and have 



A 3 (t) = -l(/8)i(t) 3 {a 3 + e(21a 2 /8 + 5a 5 )|i(t)| 2 + 0(e 2 )}. (2.52) 

This elimination of secular terms in the higher harmonics is quite general 
as also seen in Section 2.C because secular terms in the higher harmonics 
originate directly from resonant secular terms in the fundamental harmonic 
and should disappear as soon as the amplitude A of the fundamental harmonic 
is renormalized by the renormalization transformation ( |2.50| ). In this sense, we 
call secular terms in the higher harmonics as non-resonant secular terms, which 
disappear when resonant secular terms are removed by the renormalization 
transformation. 

Performing a transformation t — > t + r on ( |2.50|) and using explicit forms of 
secular terms , we obtain 

A(t + r) = A(t) + 3ea 3 \A(t)\ 2 A(t)P 1 (r)+e 2 {3al(3P 2 (r) 

+ g 1 P 1 (r)) + 10a 5 P 1 (r)}|i(t)| 4 i(t) + O(e 3 ), (2.53) 

which is an explicit representation of the Lie group underlying the renormal- 
ization transformation ( |2.50| ). Since the original equation of motion ( |2.44| ) is 

translatinally invariant , the amplitude of osillation A also enjoyes a transla- 
tional symmetry. Thus, ( 2.53| ) is also a direct result of ( 2.50| ) and the trans- 



lational symmetry. Since we could proceed this perturbation calculations to 
arbitrarily higher order although calculations becomes tedious for the higher 
order, we can also extend ( |2.53|) up to arbitrary order. Noting the important 



observation that secular terms appearing in the coefficient of the fundamen- 
tal harmonic exp(— it) in u n are polynomials of degree n and proportinal to 
|A| 2r M, we have 

oo n 

A(t + r) = A(t) + J2e n Y,Pn^\A(t)\ 2n A(t) 

n=l j=l 

= A(t) + ep 1A r\A(t)\ 2 A(t) 

+ e 2 (r 2 p 2 , 2 + rj9 2 , 1 )|i(t)| 4 i(t) + ---, 



where p n j is constant. Equating coefficients of equal powers of r in (|2.53|) and 
( [2.101 ), we have a set of equations 

— -/n\ = e n Yl e> Pn+hn \A(t)\ 2 ^A(t). (2.54) 

Some explicit forms of p n j are given as 

Pi,i = i(3/2)a 3 , P2,i = «{(15/16)a 2 + 5a 5 }, 
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p 3l i = «{(123/128)al + 5a 3 a 5 + (37/2)a 7 } 

P2,2 = "(9/8)a 2 , p 3j2 = -a 3 (45a 2 /32 + 15a 5 /2), 

p 3 ,3 = -(9i/16)a|, 

The other p nj - are obtained when higher order secular terms are calculated. 
Thus, for n — 1, eq.( 2.54 ) gives an asymptotic representation for a generator 
of the Lie group. 

= ei(3/2)a 3 |i| 2 i + e 2 i{(15/16)a 2 . + 5a 5 }|i| 4 i + ---. (2.55) 

Since Pi+jA is purely imaginary, it is easy to see that \A\ 2 is a constant of 
motion and we formally integrate ( |2.55| ) to have an asymptotic expression to 
the amplitude of the fundamental harmonic by all order: 

oo 

A(t) = Aexp[iY,£ j Pi+jM\ 2{1+j) At], ( 2 - 56 ) 

3=0 

where A is an integral constant. Then, as noted earlier, secular terms in the 
higher harmonics are removed automatically, we finally obtain an asymptotic 
solution of (|2.44|) 

u = A(t)e~ lt + eA 3 {A)e- m + e 2 A 5 {A) e - 5lt + ■■■, (2.57) 

where amplitudes of higher hamonics slave to A (e.g. see ( |2.52|) ). For the case 
that as = a-j = ■ ■ ■ = 0, the asymptotic solution (|2.57| ) may be interpreted as 
an asymptotic Fourier expansin of the Jacobi's elliptic function, which is the 
exact solution of ( |2.44j ). 

As far as the final RG equation ( j2.55|) is concerned, the following convenient 
derivation is possible without exploiting explicit forms of secular terms . Let us 
operate the linear operator L(t) defined i n (|2.30|) on the both sides of eq.( [2.53j ) 
and set r = 0. Then, taking (|2730l) and (CTT into account, we have 

L{t)A(t) = 3ea 3 |i| 2 i + e 2 {(-3/8)a 2 + 10a 5 }|^| 4 ^- (2.58) 

After (9 2 in L is eliminated iteratively , ( |2.58|) reduces to (|2.55| ). This procedure 
is valid for the secular solution ( |2.53| ) including secular terms higher order than 
0(e 2 ) in the present case. In fact, if this procedure is applied to the first order 
secular solution, we have 

LA = 3ea 3 |i| 2 i, (2.59) 

which is not correct as the first order RG equation since <9 2 in L is of 0(e 2 ) 
(see (gjjg)). 
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2.F Adiabatic Perturbation for a Soliton 



Although a soliton is a special solution of PDE, adiabatic variations of soliton- 
parameters are described by a system of ODE and we will see that a corre- 
sponding renormalization group has only one parameter. Various perturbation 
methods for solitons have been developed based on the inverse scattering trans- 
form [T3J. However, their applications are limited to nearly integrable systems. 
There are some attempts [|15|] to develope perturbation methods without utiliz- 
ing the complete integrability of the unperturbed state, which are applicable to 
a pulse solution of non-integrable systems. Here, in line with these attempts, 
we derive the well-known adiabatic equation to soliton-parameters of the K-dV 
soliton by means of the RG method. Let us consider the K-dV equation with 
a perturbation 

d t u-6ud x u + d 3 x u = ef(u). (2.60) 
Expanding u around a soliton solution 

u = s{k,z) = -2k 2 sech 2 (kz), z = x-v(k)t + £, v{k)=4k 2 , (2.61) 

we have, for the first order correction ui, 

{d t + L(z)} Ul = f{s), (2.62) 

L = -vd z - 6(sd z + d z s) + d 3 z . (2.63) 

Since k and £ in ( |2.61| ) are arbitrary parameters, zero-eigenfunctions of L(z) 
with a bounded boudary condition are given by d z s(k, z) = s', dks(k, z) = s. 
In fact, we see that zero-eigenfunctions are degenerate [16| so that 

Ls' = 0, Ls = vs', L 2 s = 0, (2.64) 

where v = dkV. Here, let us introduce an adjoint zero-eigenfunction s defined 
by the adjoint operator L to L as 

Ls = 0, s = sech 2 (kz). (2.65) 
From (|2.64j ) and (|2.65|) , we find 

POO 

< s-s' >= / ss'dz = 0. (2.66) 



DC 



In terms of these zero-eigenfunctions, we put a spatially bounded solution as 

u l = P 1 (t)8 , + Q 1 (t)s + u 1 (z). (2.67) 
Substitutiing (|2.67|) into (|2.62|) , we have an equation for ui(z) 

L(z)u! + {Pi, t + vQi)s' + Q u s = f{s), (2.68) 
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where the suffix t denotes the derivative with respect to t. Since ( 2.68|) is an 
equation for t independent and spatially bounded u\, the following consistency 
conditions must be satisfied: 

Q lt = — = constant, (2.69) 
Pi,t + vQi = constant, 

from which 

Pi,2t = -vQi,t- (2-70) 

From conditions Q2.69| ) and (|2.70| ), secular terms P\{t) and Q\{t) are found to 
be polynomials of degree 2 and 1 respectively. 

Pi = Pi,2t 2 + pi,it, Qi = q hl t, (2.71) 

where coefficients p, q are 

2pi, 2 = -vqi,i, = < *'J^J > , (2.72) 

< s ■ s > 

while pi t i is an arbitrary constant. 

Thus, a secular solution to 0(e) is found to be 

u = s(k,x + 9)+e(P 1 (t)s' + Qi(t)s), (2.73) 

where 9 — — vt + £. In order to eliminate secular terms P±,Qi by renormaliz- 
ing arbitrary soliton-parameters k and £, we introduce a renormalized soliton 
solution as 

s(k(t),x + 6(t)) = s(k, x + 6)+ e{Pi(t)s / (ife, x + 9) + Q x {t)s{k, x + 9)}, (2.74) 

where = -v(k(t)\t + f(t). Setting k(t) = k + ek^t), 9{t) = 9 + e#i(t), 
the left-hand side of (|2.74j) is expanded as 

s(k, x + 9) = s(k, x + 9) + e{9 x s' + k x s) + 0(e 2 ). (2.75) 

From ( 2.74j) and (|2.75| ), we have a renormalization transformation 

k(t) = k + eQ 1 (t;k), 9(t) = 9 + eP l (t;k) J (2.76) 

where Pi, Qi depend on k through their coefficient p, q given in ( |2.72| ). Shifting 
t to t + r in ( |2.76| ) and replacing k by k(t) — eQ 1 (t; k), we obtain a renormal- 
ization group 

k{t + r) = k{t) + eQ 1 {r;k(t)), 
9(t + r) = 9(t) - v(k(t))r + ePi(r; k(t)), 
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from which we obtain a generator of the renormalization group 



? = -v(k) + ep ljl (k), 
d 2 6 ,~ N 



(2.77) 
(2.78) 
(2.79) 



co- 



where consistency between (|2.78|) and ( |2.79|) are assured by an arbitrary ^u- 
efficient pi t \. Finally, from ( |2.72|) , (|2.77|) and ( p. 79 ), we reach the well-known 
adiabatic equations to soliton-parameters to 0(e) [nj 1 



dk 
~dt 



sech 2 (kz)f[s(k\ kz)]dx 



d 2 6 {jfy dk - dv (Ht)) 

dt 2 dt dt 



3 Reduction of Partial Differential Equations 



3. A General Discussion 



If we try to extend the Lie-group approach discussed in the preceding sections 
to PDE, there is a difficulty to calculate secular terms so that infinitely many 
divergent terms such as exponentially divergent terms as well as polynomials 
appear even in a perturbed solution of the leading order. This is due to 
the fact that the dimension of "kernels" without any boundary conditions 
of a linearized partial differential operator is infinite in general. However, 
each divergent solution belonging to "kernels" is accompanied by an integral 
constant and so unsuitable divergent solutions such as exponentially divergent 
solutions are eliminated by nullifying corresponding integral constants. Thus, 
we restrict renormalizable secular terms to polynomials even in PDE. This 
restriction is quite natural in the Lie-group appoach since a representation 
of the Lie group with multiple parameters should be expanded in terms of a 
generator as 

A(t + r,x + £) = exp(r<9t + £d x )A(t,x) 

= {l + rd t + ^d x + (rd t + ^d x ) 2 /2 + ---}A(t,x), (3.1) 
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which is a polynomial with respect to r, £ when the expansion is truncated. 
Even if secular terms are ristricted to polynomials, there are still infinitely 
many secular solutions in kernels of a linearized partial differential operator. 
For the purpose of clasifying infinitely many polynomial terms, the following 
observation is crucial. In the expansion ( |3.1|) , each monomial is accompnied 
by a diffential operator of the same order, for example, £ by d x , t 2 by d 2 
and £r by d x d t etc., while each differential operator, say, d t maps A to the 
smaller value dtA <C \A\ in the autonomous case, which will be discussed 
in the following sections. In this sense, each monomial in a representation 
of the Lie group is considered to have the same order of magnitude as the 
corresponding differential operator. On the other hand, for systems with a 
translational symmetry, each monomial with respect to t, x in secular terms in 
a perturbation series is replaced by the coresponding monomial with respect to 
r, £ in a representation of the Lie group (see, ( |2.25| )). Therefore, each monomial 
in secular terms has also the same order of magnitude as the corresponding 
differential operator. Thus, monomials in secular terms are ordered as 

x>a; 2 >s 3 , t^>tx^>tx 2 , tx > t 2 x, etc.. (3.2) 

With this kind of order for monomials in mind, we choose polynomial kernels 
from the lower order by nullifying integral constants accompanying higher- 
order kernels. In particular, the leading order secular term consists of the 
lowest order polynomial, i.e. a polynomial of degree one. As order of per- 
turbed solutions increases, the higher order polynomial kernels are taken into 
secular solutions order by order so that order of differntial operators accompa- 
nied by polynomial secular terms is kept consistent. Through this procedure, 
we can choose suitable kernels as necessary ingredients of a renormalization 
transformation. 

On the other hand, there is inevitable flexibility appearing in the choice of ini- 
tial setting of perturbation scheme, which brings about different RG equations 
from each other depending on the choice, although only one of them is often 
interesting as shown in the next example. 



3.B Nonlinear Schrodinger Equation 



As the first example, we consider a nonlinear wave equation, which is just a 
wave-equation version of the equation for a nonlinear oscillator ( 2.44J) . 

d 2 u -V 2 u + u = e m au 3 , (3.3) 

where m is a positive integer. Let us expand u around a plane-wave solution 

Wo = Aexp{i(kx — ut)} + c.c, to 2 = 1 + k 2 , (3.4) 
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as u = uq + eu\ + tui + 
For m = 1, we have 



A secular solution of 



c^iii — V 2 Ui + Mi = a«Q. 
is given by 



« i 
LP X 



3a|A|MPi(t,x,r ± )e ie - (a/8)AV ie + c.c. 
{-2iu{d t + <j0 x ) + d 2 t - V 2 }Pi = 1, 



(3.5) 

(3.6) 
(3.7) 



where 9 = kx — out, r_i_ = (y,z), Co — d^uj = k/u. As discussed in Section 
3. A, it is assumed that the leading order secular term is a polynomial of degree 
one, that is 

Pi = Pi,o,ot + Pq,i,qX + p 0>0) i ■ r±. (3.8) 
Substituting fl3.8|) into (|3.7| ), we have 

- 2za;(p lj0) o + ^Po,i,o) = 1- (3.9) 
The secular term is eliminated by means of a renormalization transformation 
A(t,x,r±) = A + 3ea|A| 2 APi(t,x,r ± ), (3.10) 
which is rewritten as, by excuting arbitrary shifts on independent variables, 

A(t + r, x + £, r ± + 77) = x, r±) + 3ea|A| 2 A(t, sc, r±)Pi(r, £, r/; £, sc, r ± ), 

(3.11) 

which gives a representation of the Lie group with multiple parameters to 0(e). 
Here, it should be noted that the coefficients of a polynomial Pi may depend 
on an "initial position", that is VjMi is replaced by pj,k,i{t, x, in Q3.11|) . It 
is seen later in ( |3.13| ), ( |3.14| ) and ( |3.15| ) that this dependence comes through 
A. In terms of a generator of the Lie group, ( |3.11| ) is expanded: 

A(t + T,x + £,r ± + 77) = exp(Td t + Zd x + Ti-V±)A(t,x,rx) 

= (l + Td t + td x + rj-V ± + ---) 

■A(t,x,r x ), (3.12) 

Equating coefficients of r, £,77 in ( 3.11 ) and ( |3.12j ), we 



where V_ 
obtain 



{d y , d z 



d t A 
d x A 
V ± A 



3ea\A\ Api A0 , 
3ea\A\ 2 Ap 0jlfi , 

o,o,i- 



3ea\A\ 2 Ap, 



(3.13) 
(3.14) 
(3.15) 



From (pHf) , ( |3~T3| ) and ( |3.14|) , Pj^i or Pj,k,i are eliminated to yield a RG equation 

(d t + ud x )A = 3iea\A\ 2 A/(2cu), (3.16) 
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which gives a steady propagation of the modulation of plane wave under effects 
of nonlinearity. 

In order to include effects of dispersion, we set m = 2 in the initial setting of 
perturbation. Then, instead of (|3.5| ), we have 



d t V - V 2 «i + ui = 0, (3.17) 

of which secular solution is 

Ul = P (t,x,r ± )e ie + c.c, (3.18) 

where Pq is the lowest order secular solution of LP Q = given by the same 
form as Pi in (|3.8| ) with a different constraint 

Pi,o,o + ^Po,i,o = 0. (3.19) 



If secular terms in (|3.18|) are removed by a similar renormalization transfor- 
mation as above, the constraint ( |3.19|) would yields a RG equation 

(d t + ud x )A = 0, (3.20) 

of which general solution is written as A = A(x', r±), x' = x — Cot. Therefore, 
it is convenient to introduce a Galilean transformation 

x = x — Cot, t = t. (3-21) 

Then, the operator L in ( |3.7| ) is rewritten and LPq = becomes 

LP = {-2iud t + (Co 2 - l)d 2 x , - Vi - 2ud t d x , + d 2 }P = 0. (3.22) 

The lowest order secular solution of ( |3.22| ) is 

Po=P { olo^' + Poli-r±, (3.23) 

where Po|i j0 and Poai are arbitrary constants. The second order secular solu- 
tion ui is given by the same form as u\ in (|3.6|) , where Pi is the second order 
polynomial solution of LPi(t,x' ,r±) = 1, that is 

n ( 2 ) 12 , (2) / , (2) 

+ P?li-r±+p?l t, (3.24) 
-2iwp$,o + 2{u 2 - l)pg j0 - 2Tr[p^ i2 ] = 1, (3.25) 

where Tr[p[) 2 g 2 ] denotes the trace of a matrix Pq\ 2 . In terms of secular terms 
( p.23| ) and ( |3.24| ), a secular solution to 0(e 2 ) is obtained as 

u = (A + eP + 3e 2 a\A\ 2 AP 1 )e w - e 2 (a/8)A 3 e 3id + c.c, (3.26) 
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which yields a renormalization transformation 

A(t, x', r±) = A + eP + e 2 3a| A\ 2 AP X . (3.27) 

Performing the shift operation t — > t + r, a;' — > x' + £, r_|_ — > r_|_ + r\ , we can 
rewrite Q3.27 ) to obtain a representation of the Lie group 

A(t + r,x' + £,r ± + r)) = A{t,x',r ± ) + eP (r,C,V,t,x',r ± ) 

+ e 2 3a\A\ 2 A{t, x', r x )Pi(r, £, 77; t, x', r ± )(3.28) 

In deriving ( |3.28| ), we have used arbitrariness of the coefficients of Pq. Again, 
note that coefficients of secular terms Pj depend on an "initial position" 

(t, x', 7*j_). As in the preceding example, this dependence comes through A 
and so the Lie group (|3.28|) also enjoys a translational symmetry which the 
original system ( |3.3j ) has. Expanding the right-hand side of ( |3.28| ) in terms of 
the generator in the same form as (|3.12 ), we have 

(3.29) 
(3.30) 
(3.31) 
(3.32) 
(3.33) 

While ( p.29|) and ( |3.30|) give order of operators d x * ~ ~ 0(e), they do not 
put any constraints on A because p|)!i,o an d Po*o,i are arbitrary coefficients of 
kernels of the linerli zed o perator L. On the other hand, ( |3.31| ), ( |3.32| ), ( p.33| ) 
and the constraint (|3.25|) yield the following nonlinear Schrodinger equation 
as a RG equation 

id t A + (u/2)d 2 x ,A + l/(2u)V\A + 3e 2 /(2to)a\A\ 2 A = 0, (3.34) 
where u = &Ioj = (1 — u 2 )/(2uj). 

It is possible to derive the final result ( |3.34| ) formally without exploiting ex- 
plicit forms of secular terms by operating the linearized operator L(r, £,77) 
defined in ( |3T22D on both sides of (|3728| ) . Noting that LP = 0,LPi = 1 
and discarding d 2 A,dtd x >A as hi gher order terms, we immediately obtain the 
nonlinear Schrodinger equation (p. 34). 



d x ,A = 


~(1) 2~(2) 
e Po,i,o 1 e Po,l,0) 


V ± A = 


f «W + e 2 r> {2) 

e /^0,0,l t t Po,0,l5 


d 2 x ,A = 


6e 2 a|i| 2 ipg j0 , 


V 2 ± A = 


6e 2 a\A\ 2 ATr[p {2 l 2 


dtA = 


3e 2 a|i| 2 i^ 2 J )0 . 



3.C Kadomtsev-Pitviashvili-Boussinesq Equation 

Let us consider a weakly dispersive nonlinear wave equation: 

d 2 u - V • {a(u)Vu - e 2 &(w)V(VV)} = 0, (3.35) 
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where the coefficient (e 2 b(u)) of a dispersive term is assumed to be as small as 
e 2 . Expanding u around constant Mo as 

u = uq + e{eu\ + e 2 u 2 + e 3 w 3 + e 4 w 4 H ), 

a(ii) = a + e[a (em + e 2 n 2 ) + e 3 {a n 3 + (a 72)M 2 } + •■■], 
b(u) = b + e6 (ewi + e 2 u 2 ) H , 

where a dash denotes the derivative with respect to u and a = a(«o) = f 2 > 
is assumed. The above setting of perturbation is justified a posteriori. We 
have the leading-order perturbed equation 

d 2 m - t> 2 VV = Lux = 0. (3.36) 

We may write a simple-wave solution of this linear wave equation in terms of 
an arbitrary function h: 

ui = h(0, (3.37) 

where ^ = x — vt. Here, it is convenient to rewrite the linear operator L by 
means of a Galilean transformation £ = x — vt, t = t as 

L(t, rr, = d 2 t - 2vd t d i - v 2 \7\, (3.38) 

where r± = (y, z), V 2 ± = d 2 + d 2 . The next order correction obeys 

Lu 2 = 0, (3.39) 

whose lowest order secular solution with respect to t, r± is given as 

u 2 (t,r ± ;0=Pol(0-r±, (3.40) 

where a monomial of t is discarded because its coefficient is constant. ^From 
the next order equation, we have 

Lu 3 = /t, / = a hh^ - bokst, (3.41) 

where the subscript £ of h, f denotes the partial derivative with respect to £. 
Since forcing terms in the right-hand side of ( 3-41|) depend only on £, ( p.41[) 
has a non-trivial secular solution with respect to t, r±: 

ft = -2vdtp§(0 - v 2 Tr\p$(a (3.42) 

where Pol(£) is a 2 x 2 marrix function of £. Summing up secular terms up to 
0(e 4 ), we obtain 

{u - u )/e 2 = h + eu 2 + e 2 u 3 . (3.43) 
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In order to remove secular terms with respect to t,r± in u 2 ,u 3 by renor- 
malizing an arbitrary function h, we introduce the following renormalization 
transformation 

h{t, r ± ; = h{i) + eu 2 (t, r ± ; + e 2 u 3 {t, r ± ; 0- (3-44) 

Since h(t, r±; £) has a translational symmetry with respect to t, rj_ as well as 
u(t, r±; £), ( |3.44| ) immediately leads to a representation of the Lie group with 
arbitrary parameters r, 77 

&(< + t, r ± + 77; = r ± ; + eu 2 (r, 77; + e 2 M 3 (r, 77; 0, (3.45) 

where it is noted again that coefficients of secular terms p depend on an "initial 
position", e.g. PoliO * s replaced by Po 2 (t, r -L> £)■ Expanding the right-hand 
side of ( |3.45| ) in terms of a generator rdt + 77 ■ Vj_ and equating the same 
monomials , we have 

Vji = epS(t,r ± ;0 + e 2 pg(t,r ± ;0, (3.46) 
V\h = e 2 Tr[p«(t,r ± ;0], (3.47) 
d t h = fptfarr,®. (3.48) 

Since Pq\ are arbitrary coefficients of kernels of the linearized operator L, 

( [3.461 ) does not put any constraints on h but provides order of the operator 
Vj. ~ 0(e). From (|3T47|) , fUS] ) and the constraint fl3gg ), we obtain the 
Kadomtsev-Pitviashvili (K-P) equation in the three dimensional space 

2vd t h£ + v 2 S/\h + e 2 (a' hh ( - b h 3 () = 0. (3.49) 

Next, we calculate higher-order corrections to the K-P equation. For u^,u^, 
we have 

Lu 4 = d 2 (a' hu 2 - b d 2 u 2 ), (3.50) 
Lu 5 = (a'o/2)(h 2 hz)z + a' Q hV\u 3 + a' \V±u 2 \ 2 

- b' (hh 3 ^-2b dlVlu 3 + g(t,r ± ;0, (3-51) 

where g(0,0;£) = 0. Explict secular solutions of ( |3.50| ) and ( p.51|) are not 
necessary for later discussions but we note 

[Lu 4 ] r±= o = 0. (3.52) 
Then, we obtain a representation of the Lie group with higher order corrections 

h(t + t, r± + 77; = h(t,rr,0 + eu 2 (r, T}; £) + e 2 w 3 (r, 77; £) 

+ e V(r,r7;0 + e 4 w 5 (r,r7;0. (3.53) 
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Since it is straightfoward but tedious to derive the K-P equation with higher 
order corrections from ( |3.53| ) by using explict forms of secular solutions, we 
follow the simplified procedure mentioned in the last paragraph in the section 
3.B. Op erati ng L(r, 77; £) on bo th sides of (|3.53| ) at (r, 77) = (0, 0) and noting 
( ED , (HI), (PI) and (pi 



we obtain the following RG equation 



(d 2 - 2vd t d^)h - V ± {(a + e 2 a' h)V ± h} 

e 2 d^{(a h + e 2 alh 2 /2)h(. - (b + e%h)h^ - 2e 2 & V 2 L /i 5 }. (3-54) 



If = ,( 3.54|) reduces to the Boussinesq equation with higher order cor- 
rections. Therefore, ( p. 54 ) may be considered as not only the K-P equation 
with higher order corrections but also the Boussinesq equation including multi- 
dimensional effects and higer order corrections. 



3.D Phase Equations 



First, we consider interface dynamics in the following genaral reaction-diffusion 
system. 

d t U = F{U) + DV 2 U, (3.55) 

where U is an n-dimensional vector and D is a n x n constant matrix. Suppose 
( |3.55| ) has an interface solution Uq = Uq(x — vt + <fi), where <ft is an arbitrary 
constant, then 

vUo, 9 + F(U ) + DU Q ,28 = 0, (3.56) 

where 9 = x — vt + <fi and the suffix 9 denotes the derivative with rspect to 9. 
Differentiating ( |3.56j ), we have 

L(9)U ,e = -(vde + F\u )-+Dd 2 )U , e = 0, (3.57) 
LU Q ,2B = F"-U* $ , (3.58) 

where F' - V — (V ■ Vu)F(U)\ u=UoJ F" ■ V 2 = (V ■ Vu) 2 F(U)\ u=Uo . It is 
convenient to introduce the following coordinate transformation 

9 = x -vt + 0, x' = x } t' = t. (3.59) 

Then, we have 

d t = dt> - vd e , d x = d x > + d 6 , 
V 2 = d 2 + 2d e d x , + V 2 x „ 

where V 2 ,/ = d 2 , + V 2 , = (0, d y , d z ). We study the solution close to U so 
that U is expanded as 

U = U + eU t + e 2 U 2 + ■ ■ ■ . (3.60) 
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The first order term obeys 

{dt> + L(9) - [V 2 X , + 2d e d x ,)D}U 1 = LU X = 0, (3.61) 
The lowest order secular solution of ( |3.61| ) is supposed to be 

U 1 = P 1 (t',r)UoM, (3-62) 
where Pi is a scalar secular function of t' and r = (x', y, z). Then, we have 

d t ,P x = 0, d x ,P 1 = 0, V^Pi = o, 

from which 

Pi = V±Pi ■ r ± , (3.63) 

where r± = (0, y, z) and Vj_Pi is a constant three-dimensional vector. 
The second order term U2 is governed by 

LU 2 = (1/2)F" • Ul (3.64) 

of which secular solution is assumed to be 

U 2 = U 2 {6) + P 2 (t', r)U ofi + R 2 (t', r)U 0>2e , (3.65) 

where U 2 is a bounded function of 9 and P 2 , R 2 are scalar secular terms. Sub- 
stituting (|3.65|) into ( |3.64|) and noting ( |3.57|) and ( 3.58|) , we obtain 



L{9)U 2 + d t/ P 2 U 0fi - Vl,P 2 DU 0>e - 2d x ,P 2 DU , 2e 
+ d f R 2 U Qj2e - V 2 x ,R 2 DU 0t2e - 2d x ,R 2 DU , 3e 
+ {R 2 - {l/2)Pt}F" ■ Ul e = 0, 

which is a equation for U 2 . Since U 2 is a function of 9 only, the following 
conditions are necessary. 

dtfP 2 = c u d x ,P 2 = c 2 , VlP 2 = c 3 , R 2 = (l/2)P? + c 4 , (3.66) 



where Cj are constant. Since Pi is a linear function of r± given in ( |3.63| ), ( p.66| ) 
yields 

d t ,R 2 = d x ,R 2 = 0, V 2 X ,R 2 = |V ± Pi| 2 . (3.67) 

Then, a bounded solution U 2 is possible only if the following compatibility 
condition is satisfied. 

< U ■ U ,9 >d t ,P 2 -<U- DU ,e > V\P 2 

<U-DU , 2e >(2d x ,P 2 + \V ± P 1 \ 2 ) = 0, (3.68) 
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where U is an adjoint function of a null eigenfunction of L and < U ■ U >= 

Now, a secular solution up to 0(e 2 ) is 

U = U (x-vt+(j ) )+eP 1 (r ± )U 0> e+e 2 {P 2 (f,x',r ± )Uo,e+R2(r±)U 0t2 e}, (3.69) 

which should be equated with the renormalized Uq(x — vt + <fi), where <p — 
4>(t', x', r±) is a renormalized phase. Setting 

4>tf,x',r ± ) = <J) + 5(t\x\r ± ), \5\ < \<j>\, (3.70) 

we expand U (x — vt + <fr) as 

U (x -vt + $) = U o (0) + 5U Qj e + (8 2 /2)U 0j2e + ■■■. (3.71) 
Equating ( 3.69|) and ( ft. 71 ), we have a renormaization transformation 



5 = ${t',x',r ± )-<f> = eP 1 {r 1 _)+e 2 P 2 {t',x',r ± ), (3.72) 
5 2 = 2e 2 P 2 (r ± ), (3.73) 

where (|3.73| ) and ( ft.66| ) are consistent with ( ft.72|) up to 0(e 2 ) if c 4 = 0. 
Following the same procedure as the previous examples, we rewrite ( ft.72j ) as 

0(t' + r, x' + Z, r± + v) = x', r ± ) + eP 1 (r t ; t\ x', r ± ) + e 2 P 2 (r, £, rj, t', x', r ± ), 

(3.74) 

from which we obtain 

= e[^(P 1 + eP 2 )] , Vl0 = e 2 [^P 2 ] o , 
d x 4 = e 2 [^P 2 ] , 8 t 4 = e 2 [d T P 2 ] , (3.75) 



where [P(t,£,t))] = P(0,0,0). Subsitution of the relations (|3775Q into ( gjjg ), 
we obtain a phase equation for the interface by replacing x' by t, x. 

(d t + 2Vd x )4> = D ± V 2 J + V\V^\ 2 , (3.76) 
D ± = <U-DU , e > I <U-U , e >, 
V = < U ■ DU , 26 > I < U ■ U 0! g > . 



If ( ft.55[ ) has a periodically oscillating solution U (—ujt + 4>) = U (9), 
the similar procedure as above leads to the isotropic Burgers equation for 
the phase 0. In this case, the linealized operator L in ( ft.61| ) is replaced by 
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L + DV 2 , L = dt — F'- and secular solutions Ui, U 2 up to 0(e 2 ) are obtained 



as 



Ui = Px(t,r)U Q (9), U 2 = P 2 (t,r)Uw + R 2 (t,r)U ,28, 
Pi = VPi-r, 

P 2 = d t P 2 t + [VP 2 ] • r + p (2) ■ rr, 
R2 = Pi/2, 

where is a 2 x 2 matrix. The constraint corresponding to ( |3.68 ) is given 
in terms of an adjoint null function U of L = d t — F'- by 

< U ■ U , e > d t P 2 - < U ■ DU ofi > V 2 P 2 - < U ■ DU^ 20 > \VPt\ 2 = 0, (3.77) 

where < U ■ U >= (U ■ U)d9 and T is a period of Uq with respect to 6. A 
representation of a renormalization group and its generator have formally the 
same form as ( |3.74|) an d ( |3.75| ), where t', Vj_ are replaced by t, V respectively. 
From ( |3.77| ) and ( 3.75| ), we obtain the three-dimensional Burgers equation 



^0 = P ± V 2 + V|V0| 2 , 
which is known as the standard phase equation. 



(3.78) 



4 Summary 



Through several examples, the perturbative renormalization group method is 
shown to be understood as the procedure to obtain an asymptotic expression 
of a generator of a renormalization transformation based on the Lie group. 
The present approach provides the following simple recipe for obtaining an 
asymptotic form of a RG equation from not only autonomous but also non- 
autonomous ODE. 

(1) Get a secular series solution of a perturbed equation by means of naive 
perturbation calculations. 

(2) Find integral constants, which are renormalized to elliminate all the secular 
terms in the perturbed solution and give a renormalization transformation. 

(3) Rewrite the renormalization transformation by excuting an arbitrary shift 
operation on the independent variable: t — *■ t + r and derive a representation 
of the Lie group underlying the reormalization transformation. 

(4) By differentiating the representation of the Lie group with respect to arbi- 
trary r, we obtain an asymptotic expression of the generator, which yields an 
asymptotic RG equation. 
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This procedure is valid for general ODE regardless of a translational symme- 
try. When the renormalization transformation is known to have a translational 
symmetry in advance, the step (3) i.e. reduction to the Lie group from the 
renormalization transformation becomes trivially simple as implied in ( 2.25Q . 
The above recipe for ODE is also applicable to autonomous PDE if we choose 
suitable polynomial kernels of the linearized operator. First, we should take 
the lowest-order polynomial, of which degree is one, as the leading order sec- 
ular term. As perturbation calculations proceed to higher order, polynomial 
kernels of higher degrees are included in the higher-order secular terms order 
by order. Through this procedure, we uniquely determine suitable polynomial 
kernels among infinite number of kernels of the linearized operator and the 
step (1) in the recipe is completed. There are no problems in the other steps. 
Thus, the present Lie-group approach is shown to be consistently applicable 
to PDE and some examples are presented. As more involved examples, we 
have succeeded in derivation of higher-order phase equations such as the non- 
isotropic Kuramoto-Sivashinsky equation and the K-dV-Burgers equation from 
the general reaction-diffusion system ( |3.55| ), which will be published elsewhere. 



Ackowledgement 



One of authers (K.N.) wishes to thank Prof.Y.Oono, University of Illinois ,for 
valuable discussions about the RG method. We also appreciate Prof. Y.Nambu, 
Nagoya University, for his introduction to the Einstein's gravitational theory. 



References 

[1] L.Y.Chen,N.Goldenfeld and Y.Oono, Phys.Rev.Lett. 73, 1311(1994). 

[2] L.Y.Chen,N.Goldenfeld and Y.Oono, Phys.Rev.E54,376(1996). 

[3] YYYamaguchi and Y.Nambu, Prog.Theor.Phys.l00,199(1998). 

[4] S. Sasa, Physica D108(1997)45(1997). 

[5] K. Matsuba and K. Nozaki, Phys.Rev.E56,R4926(1997). 

[6] T. Kunihiro,Prog.Theor.Phys.94 ,503(1995). 

[7] T. Maruo, K. Nozaki and A.Yosimori, Prog.Theor.Phys.101,243 (1999). 



31 



[8] K. Matsuba and K. Nozaki,J.Phys.Soc. Jpn.66,3315(1997). 

[9] T.Taniuti,Suppl.Prog.Theor.Phys.55,l(1977). 

[10] A. H.Nayfeh, Perturbation Met hods (John Willy and Sons, New York, 
1973). 

[11] D.V.Shirkov,intern.J.Modern Phys.A3, 1321(1988). 

[12] Y.Nambu and Y.Y.Yamaguchi, DPNU-99-09 (preprint). 

[13] for example, D.J.Kaup and A. C. Newell, Proc.Roy.Soc. London A361 
,413(1978). 

[14] for example, V.I.Karpman and E.M.Maslov,Sov.Phys.JETP 46,281 
(1977). 

[15] for example, A.Bondeson,M.Lisak and D. Anderson, Physica Scripta 
20,479(1979). 

[16] T.Ohta, Mathematical Physics in Interface Dynamics (Nihon- 
Hyouronsya,Tokyo,1997) in Japanese. 



32 



